Spontaneous-emission suppression via multiphoton 
quantum interference 



Jorg Evers^ and Christoph H KeiteP'^ 

^ Theoretische Quantendynamik, Fakultat fiir Mathematik und Physik, Universitat 

Freiburg, Hermann-Herder- Strafie 3, D-79104 Freiburg, Germany 

^ Max-Planck-Institut fiir Kernphysik, Saupfercheckweg 1, D-69117 Heidelberg, 

Germany 

E-mail: evers@physik.uni-freiburg.de, keitelSmpi-hd.mpg.de 

Abstract. The spontaneous emission is investigated for an effective atomic two- 
level system in an intense coherent field with frequency lower than the vacuum- 
induced decay width. As this additional low-frequency field is assumed to be intense, 
multiphoton processes may be induced, which can be seen as alternative transition 
pathways in addition to the simple spontaneous decay. The interplay of the various 
interfering transition pathways influences the decay dynamics of the two-level system 
and may be used to slow down the spontaneous decay considerably. We derive 
from first principles an expression for the Hamiltonian including up to three-photon 
processes. This Hamiltonian is then applied to a quantum mechanical simulation of 
the decay dynamics of the two-level system. Finally, we discuss numerical results of 
this simulation based on a rubidium atom and show that the spontaneous emission in 
this system may be suppressed substantially. 
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1. Introduction 



While spontaneous decay is a fundamental ingredient for many physical processes, many 
applications have been proposed recently where spontaneous processes are amongst 
the main limiting factors pp. These schemes usually rely on the persistence of either 
population trapped in a specific state or of coherences on timescales long as compared to 
typical atomic timescales such as the lifetime of an atomic state. A well-known example 
for this restriction is the construction of a high frequency laser, where it is hard to 
reach a population inversion as the spontaneous decay is considerable on high frequency 
transitions. For the storage and the processing of quantum information, spontaneous 
emission is a major limitation, because it is necessary to avoid all possible sources of 
decoherence in these schemes. The same holds for the secure information transmission 
using quantum effects. 
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Because of the great interest in these apphcations, various schemes to modify the 
spontaneous dynamics of an atomic system have been proposed 013111130111013 
EH HH 112] • One ansatz is the quantum Zeno effect ^ . According to the measurement 
postulate, a system is projected into one of its eigenstates upon a measurement. If the 
measurements are repeated rapidfy, the system evolution may effectively be stopped, 
because each measurement projects the system back into the initial state. The technical 
conditions on the brevity of the pulses though may not always be easily fulfilled for 
every transition, especially in the free vacuum. Another possibility is a modification of 
the vacuum surrounding the given atomic system, e.g. by an optical cavity [T. If the 
vacuum is modified such that the mode density at the frequency of an atomic transition 
is decreased, spontaneous processes on the given transition may be suppressed. Here 
the control of the environmental modes with cavities is rather challenging in reality and 
not suitable to all schemes. Thirdly it is possible to find superpositions of more than 
one upper state which are stable jl] or almost stable j3 against spontaneous decay. 
The suppression of the spontaneous decay here is due to quantum interference effects 
such as the cancellation of the amplitudes of several possible pathways between two 
system states [Zj. In spite of the conceptual beauty, the disadvantage here often is the 
difficulty to provide convenient atomic systems which fulfill all conditions such as parallel 
transition dipole moments for interference to be present. Finally, a system driven by an 
electromagnetic field periodic in time is somewhat related to systems exhibiting spacial 
periodicity such as crystals. Thus schemes relying on spacial periodicity which involve 
a change of the probability for incoherent processes such as the Borrmann effect ^Hj 
or the suppression of nuclear reactions ^3] may be transferred to laser-driven atomic 
systems jH]. Recently, this idea has also been applied to modify the decay of a three 
level system in V-configuration [H]. 

In ^U], a scheme was proposed to slow down the spontaneous decay of the upper 
state population of an effective atomic two-level system considerably. This scheme 
makes use of an intense low-frequency laser field of constant frequency and intensity 
which is applied to the two-level system such that the frequency of the low-frequency 
field is lower than the decay width of the atomic transition. This additional field 
induces multiphoton transitions between the two atomic states under consideration, 
thus allowing for alternative pathways between the two states in addition to spontaneous 
transitions. In the literature, Hamiltonians describing two-photon processes have been 
discussed to a great extend [13 HH] . However these mainly apply to systems involving 
dipole-forbidden transitions. As our aim is to inhibit spontaneous decay, it is not 
sufficient to look at the spontaneous emission on dipole-forbidden transitions as these 
rates are very low naturally. For dipole- allowed transitions, two- and four-photon 
pathways vanish due to parity reasons. Therefore it is necessary to extend the analysis 
to third-order processes ^Zj. The leading-order corrections are five-photon processes, 
which we neglect as these have a low probability in the discussed parameter range. 

Thus in this paper, we derive the Hamiltonian describing the interaction of the two- 
level system with both the vacuum modes and the additional low-frequency field mode 
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including up to three-photon processes from first principles by adiabatically removing 
the intermediate atomic states of the multiphoton processes. We calculate the coupling 
constants which were taken as free parameters in [TO] in terms of the properties of the 
atom and the laser field. Furthermore, the Hamiltonian includes Stark shift effects not 
discussed in jTO] , which however will turn out to be of no importance for the suppression 
scheme. This Hamiltonian is then applied to a quantum mechanical simulation of the 
decay dynamics of the two- level system. As it will turn out, the spontaneous decay 
may be slowed down considerably due to the additional transition pathways induced by 
the intense low-frequency field. In the final part, this is demonstrated using numerical 
results of this simulation based on rubidium atoms. 

2. Derivation of the multiphoton Hamiltonian 

The derivation starts from the usual Hamiltonian for an atomic system coupled to 
quantized electromagnetic fields ^J. The atomic system consists of two atomic states 
(typically of opposite parity) which are considered as the ground state |1) and the excited 
state 1 2) of an two-level atom, which we will denote as the effective two- level system 
throughout this paper. To model possible multiphoton transitions between the two 
effective atomic levels, unspecified auxiliary states |j) (j = 3, 4, . . .) are required, which 
will be adiabatically eliminated throughout the analysis. These states might even be 
part of the continuum; the only condition on them is that they are sufficiently far away 
from the two main atomic states as will be explained in more detail below. Since there 
are no isolated two-level atomic systems, the presence of this states is always guaranteed. 
The auxiliary states are necessary as in starting from first principles, i.e. the interaction 
Hamiltonian Hj = —dE, we only allow for one-photon transitions between the two 



Figure 1. The two- level atom and the auxiliary states |3),|4),... which will be 
adiabatically removed in the derivation of the multiphoton Hamiltonian. The auxiliary 
states are show schematically; their relative position is not restricted to be above the 
two main atomic states. After the adiabatic elimination, the effective two-level system 
exhibits multiphoton processes as symbolized in the figure. Here, the short arrows 
depict interactions with the intense low-frequency field and the long arrows interactions 
with the vacuum field. 
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effective atomic states. Multiphoton processes are then introduced by considering series 
of one-photon processes from one of the two effective states to the other effective state 
via one or more of the auxihary states, which within the adiabatic approximation occur 
simultaneously. The result will be a Hamiltonian involving the states of the effective 
two- level system only, which will turn out to be a generalization of the one used in 

The applied field consists of the low frequency field of frequency uj and the vacuum 
field modes. The following analysis is similar to the one in ^5], and the notation is close 
to [T^t im] . The system Hamiltonian is given by 

H = Ho + Hj, (1) 

Ho = nLUi\l){l\ + huj2\2){2\ + nLUj\j){j\ + nwb^b + h^Ukalak , (2) 

Hi = -dE . (3) 

Here, hui, huj2 and huj (j = 3, 4, . . .) are the energies of state |2) and auxiliary state 
respectively, is a multi index over all polarization modes and vacuum frequencies 
Uk- Throughout the analysis, we denote sums over all auxiliary states excluding the two 
effective atomic states |2) by the sum subscript j ^ {1,2}, where j is the summation 
index, b (6"'') is a low-frequency field annihilation (creation) operator, (a^) annihilates 
(creates) a vacuum photon of mode k, d is the total atomic dipole operator, and E is 
the electric field given by 



E = i [sieb - e^fet) + £k{ekak - el4) j . (4) 

The e , Cfc are polarization vectors, while the £.,£k are field amplitudes. The frequency 
corresponding to the low-frequency field operators 6, 6^ is assumed to be different 
from the relevant frequencies for the vacuum photons, which are close to the atomic 
transition frequency of the effective two-level system, such that the low-frequency 
field operators commute with the relevant vacuum field operators. In the following 
analysis, we make use of the approximation that no multiphoton transitions involving 
two or more spontaneous photons with different frequencies are included. As in our 
system spontaneous one-photon transitions are possible, transitions involving more 
than one spontaneous photon are highly suppressed; see the processes considered and 
observed in JH]- However we include multiple spontaneous photons with the same 
frequency to account for Stark shifts, which arise from interactions of the atom with the 
electromagnetical field as a sequence of an absorption and an emissions (or vice versa) 
which does not change the state of the atomic system. Also we consider all possible 
combinations of interactions with the low-frequency field as it is intense. 
The interaction part of the Hamiltonian may be expanded as follows: 

Hi = -dE = - ^ \n) {n\d\m) {m\E = - Y^dnmC^nmE 

n,m n,m 
= — {dl2(yi2 + rf21<721 j E — ^ {^ijaij + djl<Jjl + d2j<^2j + dj2(^j2^ E 
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o'nm = \n){m\ is the atomic transition operator, dnm = {n\d\m) is the dipole moment 
corresponding to transition n ^ m. Some of these dipole moments may be zero due to 
symmetry reasons, however as this depends on the configuration of the auxihary states, 
we keep this most general form of the Hamiltonian. A specific example for this will be 
given in section IH.2I where we discuss rubidium as a model system. In this expansion, 
the first part in equation (0) corresponds to direct transitions between the ground and 
the excited state by dipole coupling to one of the electric field components (i.e. one- 
photon transitions). The second addend couples both the ground and the excited state 
to one of the auxiliary states. After an adiabatic elimination of these intermediate 
states, this part will turn out to represent multiphoton transitions. Obviously n-photon 
transitions with n > 2 require transitions amongst the auxiliary states; thus these effects 
are accounted for by the last part in equation These last two parts together with 
the intense low-frequency field lead to a system evolution which qualitatively differs 
from the well-known Mollow-type evolution for two-level systems resonantly driven by 
a laser field ^B]- It is important to note that a higher order multiphoton treatment 
based on the above Hamiltonians involves operator ordering issues, as in the following 
analysis the atomic transition operators involving auxiliary states will be represented 
by expressions containing photon operators. If these issues are not taken into account, 
the resulting Hamiltonians are not self-adjoint. Thus strictly speaking, expressions like 

(dijaij + djiaji){eb - (Tb^) (6) 

in the above Hamiltonians are understood to be read as 

djiOji{eb — e*6^) + {eb — e*6^) dijOij . (7) 

As the above equation equation ^ is difficult to solve exactly, we use a perturbative 
approach which consists of an expansion in the number of interactions with the auxiliary 
levels. The lowest order Hamiltonian thus does not include any interactions with the 
auxiliary levels and may therefore be written as 

= - U,Ea,, + d,,a,,E) . (8) 



This Hamiltonian simply describes the coupling of the vacuum- and the low-frequency 
field to the effective two-level atom and does not contain contributions of the auxiliary 
states. To account for two-photon processes, we have to consider transitions from the 
ground or the excited state into and out of the auxiliary states, but no transitions 
amongst the auxiliary states. For this we introduce the transition operators <yf^-,<yf^ 
and their conjugates, where the super index (0) denotes the lowest order approximation 
whose equations of motion do not include any contribution of the auxiliary states. 
Only the operators involving the auxiliary states are approximated, as the operators 
connecting only the ground and the excited state contain the main evolution of the 
system and therefore have to be taken into account to all orders. Thus the two-photon 
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Hamiltonian is given by 



+ d2jEaf^ + dj^af^Ej . (9) 

As it will be shown in section ITT] the transition operators cr|^'*, a^2j their conjugates 
consist of sums of products of a transition operator aij {i,j G {1,2}) and a photon 
annihilation or creation operator. Thus it may easily be seen that contributions to the 
Hamiltonian in equation containing transition operators to one of the auxiliary states 
describe two-photon processes. Including up to three-photon processes, a transition 
from one of the two effective atomic levels to one of the auxiliary levels may be 
followed by a transition to another auxiliary state. Thus the three-photon Hamiltonian 
may be expressed in terms of the first order transition operators (j^j,a^j^ and their 
conjugates, whose equations of motion contain either lowest-order transition operators to 
the auxiliary state or transitions between the two effective atomic states. Also, transition 
operators between two auxiliary states are possible in lowest order aj^) ^ {1,2}), 
i.e. with no reference to the auxiliary states in the equations of motion. This yields 



i,j'^{i,2} 

- {di.Ea^l^ + djiaf^ E + d2,Ea^2j + dj^af^ E^ (10) 

i^{i,2} 

as the three-photon Hamiltonian. Here, the operator ordering was not applied to the 
term containing a double sum over j, j' for notational simplicity as it will turn out to 
be irrelevant for the present analysis. For higher order processes, 
may be obtained similarly. Introducing the coupling constants 

•1 i dij 6 i dij Cfc / 1 1 N 

Ajj = , Aijk = ^ UiJ 

where the transition dipole moments are assumed to be real, the Hamiltonians may be 
written as 

^ii-pkoton) ^ 1^^^ ^ + + ^ ^x^2k a, + A^^, al) } cri2 + h.c. , (12) 

k 

j^(2— photon) photon) 

2 

+ ^ 5Z H {\n,h + Kn^ h'^ + Y, [^-^3^ ^'^ + K^,k 4) } ^m] + h.C. , (13) 

m=l j^{l,2} k 
tt(3— photon) jjll— photon) 



m=l j^{l,2} k 

+ ^ E {^^r ^ + + E {^»'^ + 4) } 4^ + h.c. . (14) 



ij'^{i,2} 
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2.1. Atomic transition operators 

2.1.1. Equations of motion The equations of motion for the various operators involved 
may be obtained using the Heisenberg equation 

where O is an operator in the Heisenberg picture. For the transition operator a-ij with 
i = l,j ^ {1,2} - which is one of the elements that occurs in the above Hamiltonians 
in equations (fT!?| IT^ - this may be expanded as follows to the lowest two orders in the 
interaction with the auxiliary states: 



^(0) _ • , ^(0) 



1 



cr 



h 



- i uJija\/ + - [ djiau + dj2cri2 ) E , (16) 



&[f = i cuijcr^'' + ^ (^dji(Tii + dj2cri2^ ^ ^ I ^21-^4°'' - ^ dj^trf^E 1 , (17) 

\ "^{1,2} / 

where Unm = ^n—^m {n, m G N). The zeroth order equation does not include references 
to the auxiliary states, while the first order includes operators connecting one of the 
main atomic states with an auxiliary state in zeroth order. Note that for simplicity we 
omitted an addend containing the transition operator a^j in the first-order equation, as 
it will turn out to be zero (see equation (|^ ). Higher order expressions may be obtained 
similarly. For i = 2,j ^ {1, 2}, we have: 

= i ^2j(^2j + ^ (dji(^2i + dj2cr22^ E , (18) 

= i u^2iO-^j'' + ^ (dji(T2i + dj2cr22^ ~ 1 ^'^'^^^^ ~ ^ djuof^E j . (19) 

For i, j ^ {1, 2} the corresponding expression in lowest order simply becomes 

^^f = ^-^,-f. (20) 

2.1.2. First-order transition operators The basic tool to eliminate the auxiliary states 
from the equation of motion is the adiabatic integration. For this, the equations of 
motion are written in terms of slowly varying operators which are denoted by the 
corresponding symbol with a tilde and which may be seen as an interaction picture 
representation of the operator: 

^ii = o-ii, a2j = (J2j e~"^2^*, 

h = he'^\ fifc = afce'^'=*. (21) 

The adiabatic approximations with superscripts (0), (1) are transferred as their full 
counterparts. The definition of the slowly varying operators is chosen as in [13^ to allow 
for a comparison of the results. This for example yields using equations (fTT|) and 
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(HSl) 

= -i A,i an b e-i(^i^+^)* - i \*^ an e-'^"^^-^)* - i A.ifc an ak e-'("i^+"'=)* 
-i X*^i, ail al e-i('^«--^)* - i X^^ b e-'(-«+2-+'^'=)* - i X*^ b^ e-'("^^+"'=)* 

— 1 Aj2fc Cri2 Ctfc e — 1 Aj2k (^12 O-k ^ • l^^J 

Now we integrate over time, using the partial integration rule. For simplicity, only the 
first addend of equation is shown to illustrate the basic idea of adiabatic integration: 

ai°^ = - i / ^11 b e"^^"^^+^^*rft + (other addends) 

= ^^-^ ^ e-'(-^-)* - / |1 fan 1 dt 

{uJij +Uj) J {uJij +U!) [dt \ J j 

+ (other addends) . (23) 

In the adiabatic integration one now assumes that the time evolution of the slowly 
changing operators 

4 f^ii b) (24) 



dt 

(which is typically of the order of the atomic decay rate or the Rabi frequencies involved) 
is low as compared to the oscillation of the exponential function. This is fulfilled if the 
Rabi frequencies involved are not too large and if the intermediate states \j) {j ^ {1, 2}) 
are sufficiently far away from the atomic ground and excited state, which we assume in 
the following. Therefore the integral on the right hand side of equation (j^Hj) may be 
dropped, yielding using (PT|) and including all addends: 

(0) _ Ajl ^11 b ^jl ^11 b^ ,a,t , ^12 b ^-i(2<:D+a.,,)t 

■' [Ujij + UJ) [Ujij - UJ) [Ujij + ZUJ + UJk) 

^ X*2 q-12 ^_i^^t ^ Xj2k ^12 dk g-i(2c^fe+<S)t ^ \*2fc ^12 dj 

{ijJij + oJk) {uJij + OJ + 2uJk) (uJij + to) 

^ Xjik &u dk ^ \*ifc 4 gi^,t_ ^25) 

{uJij + Uk) (c^ij - l^k) 

as the lowest order adiabatic approximation for this atomic transition operator. A 
similar calculation yields 

(0) _ Aj2 ^22 , ^2 ^22 b^ Xj2k ^22 dk i^^t 

J [Uj2j + iO) [Uj2j - Uj) {Uj2j + UJk) 

I ^hk ^22 dj ^.^^^ ^ Xji a2i b ^ X]^ aai b^ ^ii^,+2^)t 

i{uJ2j - tUk) {(^2j - OJk) {uJ2j - OJk - 2u)) 

^ Xjlk 02\ dk ^iot ^ \*lfc 0^21 aj ^i(2a;fc ^26) 
{uJ2j — to) {t02j — 2uJk — Co) 

From equation (j^Uj) . we obtain for i,j ^ {1, 2} in the slowly varying operator frame: 

C=0- (27) 
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As the auxiliary states are far from resonance with the apphed frequencies, in lowest 
order of the adiabatic approximation we thus have J3] 

= = (28) 
for i,j ^ {1, 2}. Another way of seeing this is that due to the vanishing time derivative, 
in this order of approximation the populations of the auxiliary states remain constant. 
As these populations are zero if the population initially is distributed over the ground 
and the excited state, they are empty at all times. 

These results are sufficient to calculate the effective two-photon Hamiltonian 
equation Q and to obtain the first order expressions for the transition operators in 
equations (fTTjl and |T9|) . As the expressions in equations (j25|l and (j26|l only contain 
transition operators involving the ground and the excited state, the elimination of 
the auxiliary states both from the two-photon Hamiltonian and from the first order 
expressions for the transition operators is obvious. Note however that the dipole 
moments connecting the ground and the excited atomic state to the auxiliary states 
remain in the equations infiuencing the various coefficients or coupling parameters. 



2.1.3. Higher-order transition operators To calculate the first order transition 
operators, we we insert the lowest-order transition operators equations and and 
their conjugates in the equations of motion for the first order operators equations (|T7jl 
and ()19p. Then the equations are transferred to the slowly changing operator picture 
using the transformation rules equations (PT|) . The resulting expression for cr[^- , 
and their conjugates may again be adiabatically integrated yielding a^^^ a^^^ and 
a'j'2 . In addition to the corresponding operator of zeroth order, each of these operators 
contains 56 addends as first order contribution, 28 of which depend on a sum over 
auxiliary intermediate atomic states \n). For example, we have 

= ^S? + + + E (^i^(-) + • (29) 

n^{l,2} 

The operators Aij, which are proportional to cTjj, and the other first order transition 
operators may be found in Appendix A| 



2.2. Effective Hamiltonians 

In this section, we will use the transition operators derived in the last section to 
give an explicit representation of the two-photon and the three-photon Hamiltonian 
in equations and The resulting Hamiltonian will turn out to be equivalent 

to the Hamiltonian used in JU]- However here we obtain expressions for the coupling 
constants gkiQi introduced as free parameters in the Hamiltonian in ^Hj and extend 
the analysis to include Stark shift contributions to the Hamiltonian. The role of the 
Stark shifts for the decay dynamics of the effective two-level atom will be discussed in 
section 01 
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2.2.1. Two-photon Hamiltonian Using the rotating- wave approximation, we drop all 
terms oscillating with frequencies of the order of cj^ or higher in the two-photon 
Hamiltonian in equation 0. Thus we neglect counter- rotating interactions with the 
vacuum field, but do not apply the rotating-wave approximation to the interaction of 
the atom with the low-frequency field. Transferring the resulting expression back to the 
Heisenberg picture adopted in the initial equations equation (0) using the transformation 
rules in equations (EH), we obtain 

jj(2~photon) ^ ^ ^ -\- aib at (721 + 012 a,k 0-21 + h.c.j 

k 

-\-h (as al Qk + ^4 Ofc a|.) an -\- h {a^ b -\- b -\- an bb+ 6^ 6^) an 

k 

+h ^ (ttr a\ ak + as cifc al) (T22 + h^agb^ b + aw b b'^ -\- au b b -\- b^ b'^) (T22 , (30) 

k 

where the coefficients are given by ao = ^uk, 

= \^ f _ A2jAjifc _ -^2jfcA-,-i A2jfcAji \ 

\ (A - A-') A . (A - -I- ^ f A — A . — /,), -I- ,7)W ' ^ ' 



(A - Aj) Aj (Aj cjfc - cj) (A - Aj - cjfc cj) 
AjifcA2j AjifcA2j A2jfcA*i A2jfcA*i 



' Z^l (A-A,) (A, - 2cj) (A - A,- - c^fc - cj) (A, + cj^ - cu)/ ' ^ ' 



E2\jik\ljk \- 2AijfcA*ifc 
A^' 5: - (A^. + 2^,) ' ^''^ 

i^{l,2} ^ ^ ^ i^{l,2} ^ ' 

E 2Xj2kX2jk _ 2A2jfcA*2fc , , 

2A,-2A2j _ V - 2A2jA*2 



^ (A -A, + 2^) ' "-=E(a3^' (36) 

i^{l,2} ^ J ^ ^ j^{l,2} V 



«ii = y { -T^r^^^^^ - ^^^'^^'^ ) , (37) 



i^{i,2} 



(Aj + c^fc - c^) (Aj ujk + t^) 



_ / -^2j-Xj2 -^2iAj2 ^ . s 

- ^^^^ i (A^A^) + (A3A-T25) J • 

and the detunings are defined as 

A = u;2 — ^^1 — i^k + ^) ; = l^jl ~ ^k ■ (39) 

Here, the terms proportional to ao are the usual one-photon Hamiltonian parts. The 
next addends including 01,0:2 are two-photon transitions which will turn out to be of 
most interest for the following analysis. Terms with a3,...,aio may be interpreted 
as Stark shifts due to the presence of the two electromagnetic field modes. This 
Hamiltonian is an extension to the one obtained in in that it includes more 
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transitions such as direct one-photon transitions due to the vacuum and also field 
ordering effects. These ordering effects are the reason why for example and 
here are not the same other than in p[3J. The addends proportional to aii,ai2 are 
somewhat generalized Stark shifts in that they involve the population operators an and 
(J22- Usually they do not appear in effective multiphoton Hamiltonians as they may 
be dropped in a rotating wave approximation if the frequency of the involved photons 
is large enough. For example, the addend with an stems from a transition b \ 
followed by a transition b\l){j\, one of which is counter-rotating. But here b represents 
a low- frequency photon such that these terms may not be dropped a priori. However as 
discussed later the numerical simulations indicate that they do not disturb the effects 
described in ^Uj. One hint that may help to understand this is that from a quantum 
mechanical point of view, these contributions account for a distribution of the system 
state over the various Fock states of the low-frequency field without inducing atomic 
transitions. The trapping effect however does not rely on a specific configuration of the 
low-frequency field modes such as a concentration of the states to only few of the Fock 
states. 

The interpretation of the general structure of this Hamiltonian is straightforward. 
In each addend of the ai (/ = 0, ...,12), the number of A coefficients is equal to 
the number of photons exchanged. The various processes involved may also be read 
off from the A coefficients, as for example in ai in equation (p?T|l . which describes 
a transition of the effective atom from the ground state |1) to the excited state |2) 
together with the absorption of both a low-frequency and a spontaneous photon {bakO'2i 
in Hamiltonian equation ()30p). The first addend stems from a transition between |1) and 
an intermediate state \j) via a vacuum-induced transition (Xjik) and a low frequency 
field photon absorption between \j) and |2) (A2j). The third addend is due to a transition 
between |1) and an intermediate state \j) via a low frequency field photon absorption 
(Aji) and a vacuum-induced transition between |j) and |2) {\2jk)- 

2.2.2. Three-photon Hamiltonian Using the results of section !^. 1.3l in the expression for 
the three-photon Hamiltonian equation (fTIHl yields after a calculation as in section 2. II 
for the two-photon Hamiltonian the following three-photon Hamiltonian: 

^(3-photon) ^ ^i2-photon) + ^ ^ | (^^fefet + ^^fetfe + p^Jj^ + ^^^tfet) ^t^^^ + h.^. } . (40) 

k 

Thus as expected the one- and two-photon processes of this Hamiltonian are identical 
to the ones in the two-photon Hamiltonian equation (jHOI). All additional terms in 
equation (pOj) are three-photon processes. The terms with Pi, P2 are corrections to 
in the two-photon Hamiltonian, as they effectively induce the same transition. jS^ and 
are coefficients to the three-photon transition parts, which will be of most interest in 
the simulation of the decay dynamics. Third-order Stark effects vanish due to symmetry 
reasons. Introducing the detuning A„ = Uni — uj^, the explicit expressions for the P 
coefficients are given as follows (Note that the sums over j and n have been omitted 
for notational simplicity; all occurrences of these indices are summed over all auxiliary 
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Pi — . / . T — — rr + . / . 7 — ; , + 



Aj(A - Aj + O) Aj(A - Aj +u;k + 2uj) A„(Aj + Uk) 

_^ A2jA^2^j2fc ^ ^jn'^j2'^nlk ^ '^n2'^ljk'^jn 



(Aj + ojk-oo){A- Aj +0) (A - A„) (A - A^- - O) ( A„ - 20) {Aj - u) 



(A - Aj) {Aj + + 200) {Aj + UJk) {An +UJk-Uj) 

\ \* \* \ \ * \ * 



(A — Aj — UJk — 00) {Aj + 2ujk + Co") (A — Aj — uJk — ( A^- + 0;^ + 2a;) 

'^nl'^j2^jnk ^jn^j2k^nl /^-j^x 

(A - A,- - a;)(A - A„ - cjfc - cu) (A - A, - uJk){A -A^-Uk + Cj)' ^ ' 

y\jly\21'^jlk , ^n2^ij^jnk , ^jn^ljk^n2 

132 = 7-r . s.. , r + 7-r __s.. , rzrr + 



(A - A,-)(A,- + Uk) (A„ - 2a;)(A,- ^Uk- 2u) A„(A,- - u) 

_^ nl^j2k^jn ^ '^j2'^12k^2j 



(A - Aj - UJk) {A-An-ouk- 00) ( Aj - 2u!) (A - A^ + 3a;) 

+ T:J . x/. 1 TT + 



(A - A„) (A - Aj + a;) ( A^ + a;^ - 2a;) ( A„ + cj^ - a;) 
\\*\* \\*\* 

^il^21fc^il , ^i2^jnkKl 



{Aj + a;fe)(A - A^ -Uk + uj) (A - A^ + a;)(A - A„ - a;fe + a;) 

^ ^21^jlk^jl ^ ^j2^l2^2jk 



(A - Aj -ujk + uj) {Aj + 2a;fc + a;) ( A^ - 2uj) (A - A^ + a;*; + 2a;) 

+ .42) 

(Aj+a;fc-tu)(A- Aj+3a;) ' ^ ' 

\nK2Kjk . ^Ij^n2)<'jnk . \2\nK1lk 

— "TT o-\/ A + 77 r.-\/ A — ; + 



+ 



{Aj - Zui) {An - 2ui) {An - 2ui) {Aj + ujk- 2ui) (A - A„) (A - A^ - a;) 

'^2lA;iAj^^^ \2j\j2K:2k , '^12A2jA*2^. 



(A - Aj) {Aj +ujk + 2u}) {Aj - 2u}) (A - A^- + So)) ( A^ + a;^ - a;) ( A - A^ + 3a;) 



_^ \nXn\^2k ^ '^i2AnlA*„;j 



(A - Aj - UJk - 2c^)(A -An-ujk-oj) (A - A^ - a;) ( A - A„ - a;^ - uj) 



_^ ^l'j^'jnKi2k ^ ^21fc\l'^il 



{Aj -\- UJk — 2uj) {An + UJk — uj) (A — Aj — UJk — c<j) ( A^ + cj^ + 2a;) 

-^2iAjiA*ifc Ai2Aj2A2jfc 

(A- Aj -cjfe-a;)(Aj + 2a;fc + 3a;) ( Aj - 2a;) ( A - Aj + a;^ + 4a;) ' ^ ' 

r) ^12^2jfc^i2 '^12^2i'^j2fc , ^2l^jlk'^jl 

~ A / A A — ; \ + 7a — ; ^tta A — + 



Aj(A-Aj+a;fe) ( Aj + a;^ - cu) ( A - Aj + o^) (A - Aj)(Aj + a;fc) 

y^ljy^jn^n2k _|_ ^12fc^2i^j2 _|_ '^j2'^jn'^nlk 



{Aj + UJk) {An + ujk-0) Aj{A - Aj + (D) (A - A„) (A - Aj + a)) 

\*\*\* \:t:\:(:\* 

^ ^ljk^jn^„2 ^ ^2U-^jl^jl ^ ^lj^jnkK2 



An{Aj + u) {Aj + a;fe) ( A - Aj - a;^ + a)) A„( Aj + Uk) 
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Figure 2. The first four transitions involved in the coupling constant corresponding 
to the transition bbaj,ai2. The solid lines denote interactions with the low-frequency 
field, the dashed line denotes a vacuum photon emission. All paths start in state |2) 
and end in state |1). The lengths of the arrows are not related to the frequencies of 
the involved photons. 



{Aj + 2ujk -uj){A- Aj -ujk + uj) (A - Aj + uj){A - An - ujk + 00) 

(44) 



(A - A„ - cufe + uj){A - Aj -ujk + 2ij) ■ 
The structure of the coefficients is similar to the one in the two-photon Hamihonian, 
and the contributing transition pathways may be read off each addend. Figure |21 
shows the ffist four processes contributing to jS^. This coefficient corresponds to the 
transition bbalai2 in the Hamihonian equation ()4()|1 . For example, the ffist addend in 
in equation PHj) describes a path starting by a transition from the excited state |2) 
to the auxiliary state \n) induced by a low- frequency field absorption, followed by a 
transition to the auxiliary state |j) induced by a low-frequency field absorption, and 
then a vacuum photon emission with a transition to the ground state |1). The other 
transition pathways may be read off the addends accordingly. It is important to note 
that the frequency of the spontaneous photon emitted or absorbed in the multiphoton 
processes is always close to the atomic transition frequency, because the sum of the 
frequencies of all photons (with a relative sign between emitted and absorbed photons) 
is required to be close to the atomic transition frequency for a spontaneous emission to 
occur. This will be shown using a Wigner-Weisskopf-like calculation in section IHl 

For any application of the above Hamiltonian equation (jin|) . it is important to be 
able to decide which intermediate states need to be taken into account in the analysis. 
First of all, the number of exchanged photons restricts the atomic level space. For 
example in a three-photon transition from a S (angular momentum / = 0) to a P (/ = 1) 
state, states with / > 2 are not possible as intermediate states. The two important 
parameters which decide about the relevance of the remaining intermediate states are the 
transition dipole moments which connect the state in question to other relevant states, 
and the frequency separation to the other states. A reasonable parameter involving 
these quantities is e.g. the ratio Pi/ao (i = 1, . . . ,4), i.e., the coupling strength of the 
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multiphoton transitions relative to the single photon transition strength. By evaluating 
the expression for jSi for each combination of the intermediate states j and n separately, 
it is possible to compare the relative weights of the various pathways as e.g. shown in 
figure El Then the states which give rise to pathways with small relative contribution 
can be neglected. For this, it is not necessary to evaluate all possible combinations. If 
m is the principal quantum number of an intermediate state which is higher than the 
principal quantum numbers of the ground and the excited state, then the contribution 
of the corresponding intermediate state with principal quantum number m + 1 can be 
expected to be lower than the contribution of the state with quantum number m due 
to the larger frequency separation and due to a lower transition dipole moment. 



3. Decay dynamics of the effective two-level system 
3.1. General considerations 

To further understand the various contributions in the derived Hamiltonian, and to 
analyze the modification of the decay dynamics due to the additional multiphoton 
transition pathways we apply the three-photon Hamiltonian equation pn|l in a quantum 
mechanical simulation of the effective two-level atom subject to an intense low-frequency 
field and the vacuum modes. For this, we rearrange the free and the interaction part of 
the three-photon Hamiltonian as follows: 

Hq^ p/ioton) _ ^^t^ _|_ (^^flt 

k 

+ h jcji + a^b^ b + a^bb^ + ^(as «I «fc + «4 a-k al) | cm 

k 

+ h |u;2 + aQbH + aio bb^ + ^(«7 «I «fc + a/c el) | cr22 , (45) 

k 

Hi3-photon) = n J2 {(^0 + c^ib + a2b^ + PI bb^ 

k 

+p; b^ b + p; 6t 6t + p2 b fe} afc + h.c. , (46) 
^(3-photon) ^ ^ 1^^^^ bbai^+ ai2 b b ^22) + h.c. . (47) 

Hq^ photon) effective free Hamiltonian, including the Stark shifts of the two 

atomic levels. H^ac^'^"^""'^ describes the interaction with the vacuum field, i.e. the 
spontaneous decay. The contributions involving operators b and b^ in this part account 
for multiphoton transitions consisting of one interaction with the vacuum field and one 
or more interactions with the low-frequency field. Hj^^ photon) ^Qj^^g^j^^g ^g^-j^is due to the 
low frequency field alone which drive the effective system. 

In order to calculate the system evolution we notice that ^^^(^-p/ioton) pQ^^-g^j^g 
multiphoton processes only such that its coupling constants (i G {11,12}) are 
moderate even for an intense low-frequency field. Thus in the regime where the previous 
approximations such as the adiabatic elimination are valid, the evolution which gives 
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rise to the damping of the system (i.e. the couphng to the vacuum) may be evaluated 
separately from external driving fields as it is common practice in quantum optical 
calculations, see e.g. Chapter 8.6.1 in [SO]- Thus in our case, we evaluate the system 
with }j'^~P^'-°^°^) _|_ }j^^'P^°^°'^^ alone and then combine the resulting equations of motion 
with the ones resulting from the driving part photon) ^ This approximation holds 

•r W TT(-i— photon) \\ 1 1 7-7(3— p/ioton) 1 1 1 • 1 j_i j. i- ^ • jjiZ— •photon) 

if \\ril II 1 1-^0 II' which means that coupling constants m HI 

have to be low as compared to the atomic transition frequency of the effective two-level 

system. 

First, we transfer the vacuum interaction part of the Hamiltonian H^ac^'^"^"^^ to the 
interaction picture with respect to the free part H^^ photon) ^ understand the Stark 
shift contribution, we transfer the partial Hamiltonian a|, 6" (J21 as an example: 

al 6" a2i al 6" ^21 e-i*(-2-..i-....+n^) (48) 

with 

= Q-HT,i(»7+0!s)iajai-5kl) ^itJ2iio!3+Oi4) a]ai ^ (^g) 
g^ _ ^-it{ag+aio){b'<b+n) ^it(a5+a(,) (50) 

as the Stark shift contributions due to spontaneous photons Sa and due to coherent 
low-frequency photons Si,. The two exponential functions in each of the contributions 
describe the Stark shift of the upper level |2) and of the lower level |1), respectively. 
As expected, the shifts are proportional to the number of photons in the respective 
modes. The difference in the photon number factors for the two states (e.g. [b^b + n) as 
compared to b^b) is due to the fact that the example Hamiltonian part induces changes 
in the photon numbers. To evaluate their importance, these contributions have to be 
compared to the exponential factor in equation ()48j) . There is at most one photon in the 
vacuum modes {ajai G {0, 1}), thus Sa may be safely neglected. For the low-frequency 
field, it is reasonable to assume a coherent state with a large mean number of photons 
N, as this represents an intense laser field. As the relative photon number distribution 
width decreases as N~^^'^ for coherent states, and as < n < 2 in the vacuum part of 
the three-photon Hamiltonian, we approximate b^b ^ b^b + n ^ N and thus have 

SaSb ^ e-'*("«+"i«-"^-"«) ^ (51) 

as the Stark shift contribution. It is important to note that the simulated model system 
does not loose its quantum character in adopting these approximations, which are less 
inspired by physical than by numerical reasoning, as we keep the photon operators and 
the distinguishable Fock states for the electromagnetical fields. The above argument 
holds for all terms in H^ac^^°^°^^ with the same result, thus these shifts may be taken 
care of by introducing an effective atomic transition frequency 

uj = UJ2 - uJi + N {a^ + aio - «5 - "e) • (52) 

Then the vacuum interaction part becomes 

y{?-photon) = n^[ao + aib e"'^* + «2 b^ e'^* + /?* b b^ + PI b^ b 

k 

+ (3; b^ b^ e^i^* + (3lbb aai e'("-^'=)* + h.c. . (53) 
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This expression is equivalent to the interaction Hamiltonian in pj|, while we do not 
introduce the generalized ladder operators a^^^ and a^l^ {j G Z) here which were used in 
the semiclassical approximation in We proceed with the ansatz for the wavefunction 

\^{t)) = J2En{t)\2,n,0) + J2Yl Gnit)\l,n,k'). (54) 

n n k' 

Here the first index in the kets denotes the atomic state, the second slot represents 
the number of photons in the low-frequency field, and the last entry is either for the 
vacuum without photons or k for a single photon in mode k. As described earlier, we 
first derive an equation of motion for the state amplitudes due to the vacuum part of 
the Hamiltonian v}ac^'^°^^"'\ which we denote by -E^'*^ and G^''"'*'^. In a second step, we 
calculate the equations corresponding to the driving part of the Hamiltonian V^^ photon) 
with amplitudes E^. These two sets of equations of motion are then summed to give 
the equations for the full amplitudes En and G^. Inserting the ansatz equation 
in the Schrodinger equation with Hamiltonian vjac^'^"^""'^ yields as vacuum part of the 
equations of motion for the state amplitudes 



d 



k.vac 



ih jE^ = (2, n, 0|Kl'-^'^"*°"^ I*) = hJ2{i^o + PKn + 1) + f3;n) 

k 

+ (31 V(n + l)(n + 2) G'^ll' 

+ (3; ^/n{n-l) G'-^a' e^'^* } e'^^^'^'^' , (55) 
ihjGt'^^^ = {l,n,k\Vi-^'^"'"^'>\^) 

= ;i I («* + A (n + 1) + /32n) + «t K^i e'"* 
+ «; V^TTT K+i e-'"* + A Vn(n - 1) E^!^2 e''"* 

+ f33 V(n+l)(n + 2) E™^2 e'^^^*} g-'^^-^'^^* . (56) 

Formally integrating (jKH]) and inserting the result in allows to obtain an equation 
of motion for the upper state amplitudes E^"''^ of the form 

, „f n+4 

Er = - y.y. ci{t,t')Er{t') dt' (57) 



dt 



/EE c:{t,t')Erit')dt' 

"'O k l=n-4 



where Gi(t,t') are time dependent coefficients depending on the coupling constants at 
and Pj. Expanding equation (jFfj) . it turns out that all addends on the right hand side 
are of the form 

A = - / V R{uk)Xrsk ^*(a;fc)A:^,e'^^*e''^'^*'e'(^--'=+'^^)(*-*')E;;^^5(t') dt' . (58) 

"^0 k 

Here, r, s,x,y, fi,^, k and 6 are integers depending on the specific addend. R{u!k)Xrsk 
and S{ujk)\xyk are one of the coefficients ai {i = 0, 1, 2) or Pj (j = 1, 2, 3, 4) respectively, 
written as products of the dipole moment corresponding to the spontaneous transition 



Spontaneous- emission suppression via multiphoton quantum interference 



17 



Xrsk and Xxyk and of the remaining factors R{u!k) and S{ujk). In a Wigner-Weisskopf-like 
calculation, this generic contribution may be evaluated to give 



= -R{u + Kij)S*{u + ^u) p^y e^(^+'')^* - ,Jrrs{uJ + ^u)r,y{uj + ^u) ElX,{t) . (59) 

Here Tij{x) is the spontaneous decay rate of transition i j with transition frequency 
modified to x, i.e. 



where Tij is the spontaneous decay rate of transition i j with transition frequency 
LJij, and 

p^y = ^'^'^Z (60) 

is a prefactor describing the amount of quantum interference possible between transition 
r ^ s and transition x ^ y. It is zero if the dipole moments are orthogonal and reaches 
its maximum value of 1 (-1) for parallel (antiparallel) dipole moments. Thus we find 
the usual form of a square root of the product of the two corresponding spontaneous 
decay rates as characteristic for vacuum induced interference effects 0. The R and 
5* merely are prefactors which are present because the corresponding process is a 
multiphoton process. Assuming as for the Stark shift contributions that the photon 
number distribution width of the low-frequency field is negligible as compared to the 
number of photons iV, i.e. n~n-|-l ~n— 1 ~ ... ^ N, one may introduce the 
low-frequency field Rabi frequency in equation ((HZj) given by 



= 2 XijVn+l (61) 

for transition i ^ j. This finally leads to a system of equations for the upper state 
populations given by 

f^E'rit) = -Co E^{t) - ci e'^* El-\{t) - c, e"^'^* ElX.it) ' Cs e''"" El^^,(t) 

- C4 ^r+^l^^) - C5 e^^^* EZ^,{t) - c, e"^'^* El%{t) 

- C7 e^'^* El-_\{t) - C8 e-^^^* E^^,{t) , (62) 

where the Cj (i = 1 . . . 8) are constant coefficients, whose specific form which may be 
obtained by expanding equation (j57jl and (|59|l is omitted here as e.g. cq contains 
of several hundred addends in its most general form. The generalized Stark shift 
contributions in the driving part photon) Hamiltonian give rise to an equation 

of motion for the upper atomic state amplitude which is given by 

i^^Etit) = {2,n,0\vf""'''\mit)) = dse'-^'''Et_,it) + d,e-'-^'''Et^,it) (63) 

with constant coefficients di {i = 3,4). Thus the total equation of motion of the 
upper atomic state amphtude En is given by the sum of the right hand side of the 
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two contributions of equations and (jHSI) upon replacement of the amplitudes E^"-'^ 
and by 

j^E^it) = -Co - ci e'^* E„„i(t) - C2 e-'^* - (cg + it/g) e^-* E„_2(t) 

- (C4 + id,) E„+2(t) - C5 e=^'^* E„_3(t) - C6 e-=^^^* E^+s{t) 

- C7 e^'^* E„_4(t) - C8 e"^^^* E„+4(t) . (64) 

These equations are a generalization of equations (3) in ^Hj, which means that they 
have the same structure, but different coefficients Cj (z = 0, . . . , 8) and di {i = 3, 4) than 
in (TU] . The difference is due to the fact that here we include Stark shifts and do not use 
the semiclassical sum over all possible low- frequency photon numbers as in [TUj. Some 
effects of these differences will be discussed in section 13.2.21 In particular, it will turn 
out later that the driving terms which may be interpreted as generalized Stark shift 
terms do not change the decay dynamics of the simulated system in a notable manner. 

To simulate the system behavior, one has to choose an initial number of photons 
in the low-frequency mode and a possible range of deviations Ajv from this number. The 
set of state amplitudes considered in the analysis is then chosen as {En(t)\N — < 
n < N + Ajsf}. In the equations of motion for these state amplitudes according to 
equation (jMj) all references to states outside the chosen set are neglected. Thus A^r has 
to be chosen large enough such that the outermost states are barely populated during 
the calculated evolution time to avoid numerical artefacts due to the artificially added 
borders of the simulated level space. As initial condition we choose 

1 n 

|M/(0)) = |2,a,0) = -e-l-l'/2 V ^\2,n,0). (65) 

Thus the atom is in the excited state, the vacuum is assumed empty and the low 
frequency field is in a coherent state \a) which simulates a strong quantized laser field. 
The field parameter a is given by 

a = ViV e^'^ , (66) 

as |ap is the expectation value of the photon number in a coherent state. The phase ip 
accounts for the possible complexity of a. As only upper state amplitudes -E„(t) with 
A^ — AAr<?2<A^ + A7v are considered in the analysis, also the initial photon number 
distribution of the coherent state has to be restricted. Thus W is a. cutoff of the photon 
number distribution width chosen such that W < A^. Again, this avoids population 
losses at the borders of the simulated Hilbert space. P is a normalization constant such 
that (\E'(0)|\E'(0)) = 1, which is required because of the cutoff width W. As a rough 
consistency check, the range parameters A^ and W must be chosen large enough for 
the specific system parameters such that increasing these values does not affect the 
result. After solving this set of (2A7V + 1) coupled ordinary differential equations the 
total upper state population may be obtained as 

m= E (67) 

n=N-AN 
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Figure 3. Partial level scheme of rubidium taken as the simulated model system. 
The spontaneous emission of transition |2) ^ |1) is to be slowed down by the discussed 
scheme. The figure is not drawn to scale. 



This population may then be compared to the exponential decay for the two-state system 
without an additional low-frequency field. 

3.2. Rubidium as an example 

3.2.1. Model system To demonstrate the feasibility of the scheme, we use rubidium as 
our model system. The simulated atomic levels are shown in figure 01 The population 
is assumed to be in the 5P3/2 state initially. Without additional fields, this state decays 
to the 55'i/2 state with a decay rate of = 37.5 ■ 10^ PT|. As discussed in the 
previous section, within the adiabatic approximation the intermediate states are never 
populated, see equation (j2HI)- Thus it is possible to include states even if there are decay 
channels leading out of the simulated level space. Also, low spontaneous decay rates 
do not in general mean that the corresponding transitions do not need to be taken into 
account, as the low rates may be due to the small energy spacing between two states. 
The dipole moments however which are important for the low-frequency field may be 
comparable to the other transition dipole moments. In the simulation, only dipole- 
allowed transitions are considered. This amounts to a simplification of the equation of 
the upper state amplitudes equation (j64|l by e.g. eliminating ci and C2 which rely on 
dipole-forbidden transitions, as in addition to the spontaneous photon a second photon 
is exchanged with the low-frequency field. The energies of the various states are given 
by 

^1 = eV , 82 = 1.589 eV , 

£3 = 2.496 eV , ^4 = 1.560 eV , 
85 = 2.400 eV , 86 = 2.400 eV , 
87 = 2.940 eV , 87 = 2.950 eV . 

The energy separation {85 — 8q) is about 55 MHz; the main transition has a frequency 
of UJ21 = 3.84 ■ lO^'' Hz. We consider decay rates as obtained by theoretical calculations 
in [21] , where also the required branching ratios of the various decay pathways from the 
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Figure 4. Dependence of the total upper state population on fli2. The chosen values 
are fli2 = 1; 3; 5; 8; 10; 12 [lO^^Hz]. The graph order corresponds to the Rabi frequency 
where the highest graph corresponds to the largest driving strength. The graph for 
r2i2 = 1 • 10^^ Hz is almost on top of the dashed reference curve. 
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For the transitions which are hsted in |22I; the above values are in reasonable agreement 
with the data reported there. The system parameters for the figures are chosen as 
N = 10^ Ajv = 15000, W = 500, = 0, = 0.1 MHz and p^^' = 1 if not stated 
otherwise. Here N only aff'ects the initial population of the low-frequency Fock states. 
It is chosen independent of the field strength of the low-frequency field in order to 
estimate the dependence of the population dynamics on the initial conditions. The value 
p^'f = 1 has been chosen for simplicity, but is not required to for our scheme to work. 
A calculation with a p-value of 0.5 reduces the trapping duration e.g. of the top curve 
in figure m by less than a factor of 3, such that the remaining effect is still considerable. 
For this one should note that other than in most previously studied systems exhibiting 
quantum interference effects |11 El El IZ|, in our setup non-zero values for can be 
found in any atomic system, which is due to the fact that here the two transitions do 
not need to have a similar transition frequency or a common atomic state. To speed 
up the numerical calculations, the driving terms in equation are suppressed by the 
replacement 



excited states are given: 



di^ gdi{i = 3, 4) 



(68) 
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Figure 5. Dependence of the population trapping on the low-frequency field frequency 
u!. Chosen values are w = 10^; lO'^; 10^; 10^; 10^ 10^ Hz. (a) Co = f 0^; IQ-* Hz; (b) 
w = f 0^ Hz ; (c) w = f 0"^ Hz and (d) w = 10^; 10^ Hz. The two curves in (d) can be 
distinguished as the one corresponding to the lower frequency uj slightly wiggles. The 
dashed curve is the reference. 

in equation ()64|) . where we choose g = 1/1000 in the following calculations. For this one 
should note that for the above parameters the simulation consists of 2Ajv + 1 = 30001 
coupled complex differential equations. The validity of this suppression is discussed in 
section 13.2.21 In the following figures EHZl dashed lines are reference curves given by 
exp(-f2it). 

3.2.2. Results The first figure El shows the total upper state population n(t) for 
different values of Qi2. The chosen values for this Rabi frequency are Q12 = 
1; 3; 5; 8; 10; 12 [10^^ Hz]. As expected an increasing Rabi frequency of the low-frequency 
field increases the amount of trapping. The reason for this is that the relative probability 
of low-frequency-field-assisted transitions increases with an increasing field strength. 
However the Rabi frequency must not be chosen too high as otherwise the adiabatic 
approximation is not valid any longer. Also higher-order processes with more than 
three exchanged photons have to be considered if the Rabi-frequency is chosen too 
large. The Rabi frequency Q12 = 10^^ Hz chosen in most of the figures is about 1/40 
of the atomic transition frequency 0021, which is well within the validity range of the 
applied approximations. 

Figure El shows the role of the low-frequency field frequency Co. The chosen values 
are ui = 10^; 10^; 10^; 10^; 10^; lO^Hz, while the other parameters are chosen as in figurelH 
with Q12 = 10^^ Hz. As long as u is low as compared to the natural decay width ri2, 
the result is independent of u and the population plots are on top of each other. For 
higher frequencies, the population trapping decreases until the oscillation due to u is 
visible for u ^ ri2. The behavior of the system in the limit of small frequencies u is 
somewhat different from the behavior found in ^U]. There the trapping was found to 
improve with decreasing u until in the limit a; — the decay was completely stopped. 
As already discussed in |in|, this is the expected behavior for the system simulated in 
pn] which is equivalent to a system with near-degenerate upper states (SBj ■ But while in 
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Figure 6. Dependence of the population trapping on the phase of the initial low- 
frequency field coherent state. The chosen values are = —0.25; 0; 0.25; 0.5; 0.75; l[7r]. 
Plots are (a) for (f) = 0.5tt and the dashed reference, (b) for cf) = — 0.257r; 0.757r, (c) for 
(j) = 0.257r, and (d) for = 0; tt. 

our model in jTUj all other parameters were kept constant in changing u, in the present 
analysis other parameters such as the various coupling strengths also change with ui, 
thus leading to a different behavior. Still the slowing down of the spontaneous emission 
is most pronounced for low values of u. However this is a key ingredient of all similar 
quantum interference effects j23] which may be depicted as follows: The smaller the field 
frequency is, the harder is it to distinguish between the various interference pathways, 
which leads to stronger quantum interference effects. As the scheme relies on the fact 
that photons with nonzero frequency may be exchanged during atomic transitions, the 
singular case of zero frequency, i.e. a static field, is excluded from our analysis as in 

m- 

In figure El the phase of the initial low-frequency field coherent state in 
equation is varied. The chosen values are (p = —0.25; 0; 0.25; 0.5; 0.75; l[7r]; the 
other parameters are chosen as in figure ID with Q12 = 10^^ Hz. The phase of the 
initial coherent state crucially influences the decay dynamics of the effective two-level 
system. By choosing the phase, the effective spontaneous decay varies between the usual 
decay rate and the maximum trapping for (p = 0; tt. This may be understood from 
equation (j6l|l which shows that the decay dynamics of the state amplitude En depends 
on the neighboring amplitudes En±m {m = 1, . . . ,4) which change relative to En with 
varying phase 0. 

Figure U\ (a) shows that the initial photon distribution width of the coherent low- 
frequency laser field does not influence the result of the numerical calculation notably. 
This is a consistency check of the approximate coherent state in the ansatz for the 
wavefunction equation (jHS). The chosen values are W = 100, 300, 500, 1000, 2000, 3000; 
the other parameters are chosen as in figure lU with Q12 = 10^^ Hz. 

Figure U\ (b) shows that the number of photons of the coherent low-frequency laser 
field does not visibly influence the result of the numerical simulation. To evaluate the 
influence of the initial conditions, this value was chosen independent from the intensity of 
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Figure 7. Consistency checks of the simulation, (a) Dependence of the population 
trapping on the initial laser distribution width W. The chosen values are W = 
100,300,500,1000,2000,3000. (b) Dependence of the population trapping on the 
number of photons in the low-frequency mode considered in the simulation. The 
chosen values are N = 10^, 10^ 10^ 10^ 10*, 10^. (c) Dependence of the population 
trapping on the driving strength considered in the simulation. The chosen values 
are g = 10~^; 10""*; 10~^; 10~^; 10~^; 1. (d) Dependence of the population trapping 
on the number of simulated photon number states. The chosen values are Ajv = 
3000,5000,8000,10000,12000,15000. In all subfigures, all plots but the dashed 
reference are on top of each other. 

the low-frequency field. The chosen values are = 10^, 10^, 10^, 10^, 10^, 10^; the other 
parameters are chosen as in figure H] with Q12 = 10^^ Hz. Increasing broadens the 
initial photon number distribution. However if is large enough, neighboring photon 
number states have similar amplitudes up to a possible phase. This does not change 
notably for increasing A^, and as the equation of motion for state amphtude En depends 
on the neighboring amplitudes En±m (m = 1, . . . , 4) as discussed before, the results are 
independent of A^ if it is not too small. 

Figure [7| (c) shows the dependence of the population trapping on the damping 
factor g of the driving terms which was introduced in equation (jHEI) to speed up the 
numerical calculations. The chosen values are g = 10~^; 10~^; 10~^; 10"^; 10~^; 1; the 
other parameters are chosen as in figure lU with = 10^^ Hz. The Figure does not 
show a dependence on the damping of the driving strength. This may be explained 
along the lines of the interpretation of figures [7| (a) and (b), as these driving terms 
merely account for a broadening of the photon number distribution of the low-frequency 
field. The trapping mechanism however does not rely on a specific distribution width. 
Also, the driving terms are due to the same interactions as the ones which lead to 
the Stark shift discussed before. The last consistency check is shown in figure [7| (d), 
where the number of simulated photon number states A^v is varied. The chosen values 
are A^v = 3000,5000,8000,10000,12000,15000; the other parameters are chosen as in 
figure m with Q12 = 10^'^ Hz. Again, there is no visible dependence on A/^r for the chosen 
values which shows that A^v = 15000 is large enough to eliminate errors due to border 
losses. 
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4. Summary 

In summing up, we have derived an explicit expression for the multiphoton Hamiltonian 
describing the interaction of an atomic two-level system with both the vacuum field and 
an additional intense low-frequency laser field including up to three-photon processes. 
This Hamiltonian was used in a quantum mechanical simulation of the decay dynamics 
of a two level atom subject to the intense low- frequency laser field and the vacuum field 
modes. Using this simulation it was shown that the usual spontaneous decay found on 
one of the transitions in rubidium may be decelerated considerably by a suitably chosen 
low-frequency field. 
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Appendix A. First order transition operators 

In this appendix explicit expressions for the first-order transition operators calculated 
in section 12.1.31 are given. These transition operators may be written as 

n^{l,2} 
n^{l,2} 

The corresponding operators crj}'' and may be obtained by conjugation. The 
operators A^j {i = 1, 2) are proportional to aki- They are given by 

A\]in)=Ail,n), (Al) 

where 

J^^Q, _ '^jnk^nak "^fc ^k <^ aa _^ ^nak^'jnk ^k ^fc (^aa 

[U^ - UJj){uJc, -UJn- ^k) iyJa " ^j){!^a - + ^k) 



iOn — UJj){uJa — UJn — Uj) {iOa — UJj — lijj) {uj„ — UJr, — 00 ] 
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